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FOREWORD 


The  work  presented  in  this  report  was  performed  by  Kaman 
Aerospace  Corporation  under  Contract  DAAJ02-70-C-0012 
(Task  1F162204AA4301)  for  the  Eustis  Directorate,  U.  S.  Army 
Air  Mobility  Research  and  Development  Laboratory,  Fort 
Eustis,  Virginia.  The  program  was  implemented  under  the 
technical  direction  of  Mr.  Joseph  H.  McGarvey  of  the 
Reliability  and  Maintainability  Di-vision*  and  Mr.  Arthur  J. 
Gustafson  of  the  Structures  Division.**  The  report  is  pre¬ 
sented  in  four  volumes,  each  describing  a  separate  phase  of 
the  basic  theory  of  structural  dynamic  testing  using  impedance 
techniques . 

Volume  I  presents  the  results  of  an  analytical  and  numerical 
investigation  of  the  practicality  of  system  identification 
using  fewer  measurement  points  than  there  are  degrees  of 
freedom.  The  parameters  in  Lagrange's  equations  of  motion, 
mass,  stiffness,  and  damping  for  a  mathematical  model  having 
fewer  degrees  of  freedom  than  the  linear  elastic  structure  it 
represents  may  be  determined  directly  from  measured  mobility 
data.  Volume  II  describes  the  method  of  system  identification 
wherein  the  necessary  impedance  data  are  experimentally  deter¬ 
mined  by  applying  a  force  excitation  at  a  single  point  on 
the  structure.  Volume  III  presents  a  method  cf  determining 
the  free-body  dynamic  responses  from  data  obtained  on  a  con¬ 
strained  structure.  Volume  IV  describes  a  method  of 
obtaining  the  equations  for  the  combination  of  measured 
mobility  matrices  of  a  helicopter  and  its  subsystems.  The 
response  of  the  combination  of  a  helicopter  and  its  sub¬ 
systems  is  determined  from  data  based  on  the  experimental 
results  of  the  main  system  and  subsystems  separately. 


^Division  name  changed  to  Military  Operations  Technology  Division. 
**Division  name  changed  to  Technology  Applications  Division. 
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INTRODUCTION 


The  dynamic  testing  of  helicopters,  both  full  scale  and 
model,  is  required  to  anticipate  the  response  of  the  vehicle 
to  the  vibratory  load  spectrum  to  which  it  is  subjected  in 
flight.  These  responses  were  computed  from  an  analytical 
model  which  is  partly  developed  from,  and  verified  by,  data 
obtained  from  testing,  such  as  resonant  frequencies,  res¬ 
onant  responses,  and  damping. 

The  uncertainties  in  the  experimental  determination  of  mode 
shapes  and  resonant  frequencies  for  unconstrained  structures 
relate  to  the  inherent  difficulties  involved  in  supporting 
the  structure  so  as  to  simulate  free-body  conditions. 

During  vibration  tests  the  model  or  full-scale  vehicle  must 
be  supported  in  a  manner  to  prevent  interference  with  the 
response  of  the  structure.  In  order  to  simulate  the  free- 
body  boundary  conditions  of  a  helicopter  in  flight,  the 
usual  procedure  has  been  to  support  the  structure  on  a 
system  which  is  relatively  soft  to  negate  the  effect  of 
rigid-body  modes  on  the  elastic  modes  of  the  helicopter. 

A  commonly  used  technique  consists  of  supporting  the  fuselage 
on  cables  and  bungee  attached  to  the  hub.  As  helicopters 
become  larger  and  heavier,  it  will  be  necessary  to  construct 
massive  structures  capable  of  supporting  the  total  weight  of 
the  vehicle.  This  situation,  which  is  expensive  and  requires 
substantial  design  effort,  will  continue  to  deteriorate. 

In  addition,  there  is  some  uncertainty  in  the  effects  of  the 
suspension  dvnamics  and  nonlinearities  on  the  helicopter 
response . 

A  practical  method  which  eliminates  the  need  for  soft  sus¬ 
pensions  in  order  to  simulate  in-flight  boundary  conditions 
will  be  of  great  help  in  the  dynamic  testing  of  helicopters. 
Such  a  method  would  yield  better  correlation  between  experi¬ 
mental  end  flight  results  and  significantly  reduce  the  cost 
of  this  testing. 

The  method  described  in  this  report  has  several  desirable 
features.  It  uses  the  measured  forces  of  constraint  to 
convert  the  measured  constrained  helicopter  responses  to 
free-body  responses.  The  structure  being  tested  is  con¬ 
sidered  to  be  supported  on  real  supports,  but  the  character¬ 
istics  of  the  supports  themselves  are  not  required  since 
only  their  measured  reactions  are  used.  The  procedure  uses 
only  data  which  is  actually  measured,  and  no  quantitative 
assumptions  are  employed. 
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THEORY 


For  sinusoidally  varying  forces  at  a  frequency  w,  applied  at 
M  points  on  a  structure,  a  vector  (?\  }  can  be  defined  which 

represents  the  complex  amplitude  of  applied  force  at  each 
of  the  points.  Similarly,  (y ^ }  is  a  vector  representing  the 

conplex  amplitudes  of  the  deflection  at  each  of  N  points 
resulting  from  the  force  vector  (fj).  There  is  no  necessity 

for  the  force  excitation  points  to  coincide  with  the  points 
at  which  the  deflection  response  is  measured.  The  complex 
displacement  mobility  matrix,  [Y] ,  of  order  NxM  represents 
the  relationship  between  the  applied  forces  and  responses. 
The  vectors  { f ^ } ,  (y^ >  and  the  matrix  [Y]  and  u,e  relation¬ 


ship  among  them  are  written 


{y  }  =  1 
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The  previous  expressions  are  not  restricted  to  displacement 
response;  exactly  the  same  relationships  apply  for  velocity 
and  acceleration  response.  The  displacements  can  also,  with 
no  modification  in  the  analysis,  represent  displacements  or 
rotations  in  two  or  three  directions  at  one  geometrical  point 
by  alloting  one  element  in  each  vector  for  each  generalized 
displacement.  Similar  considerations  apply  to  the  forces 
or  moments.  There  is  no  necessity  for  the  matrix,  [Y]  ,  to  be 
square;  ic  will  contain  one  row  for  each  measured 
displacement,  one  column  for  each  po:'nt  at  which  an  excitation 
is  applied. 


Consider  a  matrix,  [F]  ,  containing  several  applied  vectors 
and  a  matrix,  [y] ,  containing  the  corresponding  deflections, 
as  follows: 


?*’]  =  [7, 
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!  2: 


•  •  •  •  ] 


[y]  =  iv1  | 


•  • 


.] 


2 


V 


and .then 


[y]  =  [Y][F] 


(2) 


If  [F]  is  a  nonsingular  matrix,  then  the , desired  result,  the 
response  of  the  points  of  interest  to  single  forces,  may  bo 
written 

[Y]  =  [y] [F]"1  (3) 

In  Equation  (3) ,  the  matrices  [y ]  and  [F]  consist  of  measured 
data.  When  the  externally  applied  loads  only  are  included 
in  [F] ,  then  [Y]  is  the  mobility  of  the  structure  as  tested 
on  the  actual  supports.  If  [F]  includes  any  of  the  forces  of 
constraint,  then  [Y]  is  the  mobility  of  the  structure  with 
those  constraints  removed.  If  all  the  forces  of  constraint 
are  included  in  [F] ,  then  [Y]  is  the  mobility  of  the  uncon¬ 
strained  system  or  free  body. 

As  indicated  previously,  [F]  must  be  nonsingular  and  thus 
have  an  inverse.  If  there  are  M  forces  to  be  considered, 
including  the  forces  at  the_constraints ,  then  M  sets  of 
vectors  of  applied  force,  { f ^ } ,  must  be  applied  and  all  of 

these  vectors  must  be  linearly  independent.  The  independence 
criterion  may  be  acnieved  by  applying  an  external  force  at 
each  constraint. 


FORCES  AT  CONSTRAINTS 


If  an  exciting  force  is  applied  at  the  j-th  constraint,  the 
force  vector  will  be  of  the  form 


<V  - 


f  +  f 
rj 


(4) 


where  the  f  are  the  forces  at  the  constraints  due  to  the 

rj 

applied  force  f. 


For  an  optimally  designed  support  system,  the  maximum  number 
of  constraint  forces  is  six  and  the  [F]  matrix  relating  the 
support  forces  will  be  of  the  following  form: 


( 


! 
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where  the  elements  in  each  column  represent  the  actual  meas¬ 
urements  taken  for  the  force  applied  at  each  of  the  constraints. 
Because  the  loads  are  independent,  the  preceding  matrix  will 
always  be  nonsingular. 

In  conjunction  with  the  force  measurements,  the  displacements 
(or  accelerations)  are  recorded  at  the  points  of  interest  on 
the  structure  and  one  column  of  Ey ]  is  formed  for  each 
column  of  [Fc] . 

Then,  as  in  Equation  (3)  , 

[YFB]  =  [y][Fc]_1  (6) 


where  [Ypfi]  represents  the  deflection  of  each  point  of  in¬ 
terest  on  the  structure  due  to  each  of  the  applied  loads  at 
the  supports,  and  is  the  free-body  mobility  matrix.  The  [F  ] 
matrix  was  obtained  at  one  particular  frequency,  and  the 
same  procedure  must  be  implemented  over  the  frequency  band¬ 
width  of  interest. 

Normally  it  is  required  to  determine  the  response  of  the 
helicopter  due  to  force  excitation  at  points  other  than  the 
supports.  Therefore,  the  structure  must  also  be  forced  at 
these  points  and  the  forces  at  the  constraints  must  be 
measured.  If  [Fc  ]  is  a  matrix  representing  the  forces  of 

constraint  for  each  point  of  excitation  and  [FAJ  is  a  matrix 

representing  the  applied  loads,  then  the  [F]  matrix  becomes 


However,  if  unit  forces  are  applied,  the  matrix  [F  ]  becomes 
the  unit  matrix  and  Equation  (7)  becomes 
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(8) 


[F]  = 


[i]  !  o  1 


The  inverse  of  this  matrix  is  readily  obtained  and  is  given 
by 


[F] 


-1 


■t 


.in. 


ltPc]'1[PcAl  1  [Fc! 


"I 


(9) 


or  by  simply  inverting  the  [F]  matrix  in  its  entirety  by 
numerical  means. 

ERROR  ANALYSIS 

The  technique  of  structural  dynamics  testing  utilizing 
impedance  methods,  to  be  of  any  practical  engineering 
significance,  must  be  functional  with  a  reasonable  degree 
of  experimental  error.  Measurements  of  the  complex  ampli¬ 
tudes  of  displacement  required  to  implement  the  method  will 
be  subject  to  experimental  errors  of  various  types.  In 
general,  all  errors  can  be  classified  as  either  random  or 
bias. 


In  the  present  analysis ,  the  simulated  test  data  were 
polluted  with  measurement  error,  and  computer  experiments 
were  conducted  to  test  the  sensitivity  of  the  method  to 
error.  The  20-degree-of-f reedom  "beam-type"  representation 
of  a  helicopter  and  the  1 8-degree-of-f reedom,  three- 
dimensional  mathematical  model  were  analyzed  incorporating 
measurement  error.  The  simulated  test  data  in  the  form  of 
acceleration  mobility  and  reaction  force  measurements  were 
both  polluted  with  a  5  percent  bias  error  and  a  +5  percent 
random  error  on  amplitude  and  a  +1  degree  error  on  phase 
angle.  The  random  amplitude  phase  angle  errors  were  assumed 
to  be  uniformly  distributed.  The  resulting  distribution  of 
random  error  between  the  selected  limits  is  conservative 
compared  to  limits  applied  to  a  normal  distribution. 
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V 


APPLICATION  OF  THE  THEORY 


The  method  has  several  attributes  which  make  it  especially 
suitable  for  practical  application,  including  the  use  of  only 
measured  data,  the  lack  of  quantitative  assumptions  regarding 
the  mass  or  stiffness  distributions  or  the  assumption  of 
infinitely  rigid  supports.  There  are,  however,  as  in  all 
procedures,  certain  considerations  involved  in  planning  an 
efficient  and  accurate  application  of  the  method. 


At  each  frequency  it  is  necessary  to  conduct  one  test  for 
each  constraint.  It  is  possible  to  constrain  all  rigid- 
body  motions  with  six  constraining  forces.  There  is  no 
necessity  for  such  complete  constraint,  however.  During 
the  Resign  of  a  test,  test  configurations  should  be  considered 
which  tr’JLow  freedom  of  motion  in  the  horizontal  plane  and 
around  t he**^ rt i ca i  axis.  Such  an  arrangement  for  a  typical 
helicopter  migi>*  consist  of  three  supports  at  the  base,  each 
resisting  only  ve^i  *al  motion  and  mounted  on  a  device  which 
allows  unconstrained,  Jjj*  limited,  motion  in  the  horizontal 
plane.  For  sinusoidal  force  ^citation  at  the  rotor  hub, 
it  is  necessary  to  measure  the®^mplex  amplitude  of  the 
response  of  the  points  of  interest  on  the  structure  as  well 
as  the  forces  induced  at  the  constraint  points.  It  is  also 
necessary  to  apply  vertical  excitation  at  each  support  and 
measure  each  of  the  vertical  forces  of  constraint.  In 
addition,  any  other  shaker  position  or  orientation  could  be 
used  while  the  three  vertical  forces  are  measured. 


Theoretically,  the  stiffness  characteristics  of  the  supports 
are  immaterial.  However,  these  characteristics  do  affect 
the  magnitudes  of  the  forces  and  responses  which  will  be 
measured.  The  performance  of  transducers  such  as  accelerom¬ 
eters  or  load  cells  is  dependent  on  the  magnitude  and 
frequency  of  the  quantity  being  measured.  Thus,  for  the 
most  reliable  results,  the  supports  should  be  designed  and 
the  transducers  selected  for  optimum  operation. 

A  computer  simulation  of  alternative  test  configurations 
can  be  an  invaluable  aid  in  the  preliminary  design  stage  of 
an  actual  test  program.  Prior  to  implementation  of  a  test 
plan,  a  computer  simulation  applied  to  this  method  using  an 
intuitive  analytical  model  of  the  structure  can  be  used  to 
determine  the  sensitivity  to  error  and  the  expected  accuracy 
of  the  results  of  the  various  arrangements  considered.  It 
should  include  realistic  experimental  errors,  approximate 
constraint  characteristics  and  a  frequency  spectrum  bounding 
the  range  of  interest. 
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COMPUTER  SIMULATION  RESULTS 


i 


A  computer  simulation  study  was  conducted  to  determine  the 
practicality  of  the  theory  described  in  the  subject  report. 

The  digital  computer  program  listing  in  FORTRAN  IV  language 
and  a  description  of  the  program  input  cards  are  presented  in 
the  appendix.  Two  basic  configurations  were  analyzed,  a  20- 
degree-of-freedom  "beam-type"  representation  of  a  helicopter 
and  an  18-degree-of-f reedom,  three-dimensional  structure. 

Table  I  presents  a  description  of  the  "beam-type"  model 
including  mass  and  stiffness  distributions.  The  discrete 
masses ,  coordinate  locations ,  member  internal  spring  rates 
and  ground  springs  are  given  in  Table  II  for  the  18-degree- 
of-f  reedom  analytical  model.  For  the  models  investigated, 
the  simulated  test  data,  acceleration  mobility  amplitude 
and  reaction  forces,  were  separately  polluted  with  measurement 
errors  of  5  percent  bias  error,  +5  percent  random  error  and 
+1  degree  error  on  phase  angle.  Tables  III  and  IV  present  a 
compilation  of  a  portion  of  the  computer  experiments  performed 
in  implementing  the  theory  presented  in  this  report. 

The  effect  of  measurement  error  on  the  free-body  real  and 
imaginary  acceleration  frequency  response  for  the  20-degree- 
of-f reedom  "beam-type"  representation  of  a  helicopter  is 
shown  in  Figures  1  and  2,  respectively.  In  each  instance 
the  exact  curve  represents  the  result  of  the  analysis 
utilizing  simulated  experimental  data  with  zero  error. 

Figures  3  and  4  present  similar  data  for  the  18-degree-of- 
f reedom,  three-dimensional  mathematical  model.  As  evidenced 
by  the  aforementioned  figures,  the  theoretical  development 
presented  herein  is  essentially  insensitive  to  the  dearee  of 
error  incorporated  in  the  analysis. 

Figures  5  and  6  portray  the  results  of  implementation  of  the 
method  using  the  measured  forces  of  constraint  to  yield  the 
free-body  response  of  the  20-degree-of-f reedom  "beam-type" 
structure.  Two  simulated  constrained  modes  are  shown  on  the 
figures  with  the  specimen  supported  on  real  supports  of 
10,000  lb/in.  and  50,000  lb/in.  spring  rates,  respectively. 

The  free-body  response  of  the  structure  achieved  by  deleting 
the  effects  of  supports  for  both  constrained  boundary  con¬ 
ditions  is  also  presented.  It  is  important  to  note  that 
regardless  of  support  system,  application  of  the  method  yields 
identical  unconstrained  or  free-body  boundary  conditions  for 
the  structure.  Many  computer  experiments  were  conducted, 
including  the  effect  of  measurement  error,  to  obtain  the 
data  presented  in  Figures  5  and  6.  Comparison  of  the  simu¬ 
lated  "exact"  free  response  and  the  response  obtained  incor¬ 
porating  realistic  measurement  error  indicates  the  results 
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TABLE  I.  2 O-DEGREE-OF- FREEDOM  MODEL,  "BEAM  TYPE 


TABLE  III .  20-DEGREE -OF- FREEDOM  11  BEAM  TYPE”  CONFIGURATION 


**  Applied  forces  at  Sta  1,  5,  10  and  18 


TABLE  IV.  THREE-DIMENSIONAL  CONFIGURATION 


FREE-BODY  ACCELERATION  MOBILITY,  REAL  (G/LB-FORCE) 


Figure  1.  The  Effect  of  Error  on  "Beam-Type"  Model 
Driving  Point  Real  Acceleration  Response 
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FREE-BODY  ACCELERATION  MOBILITY,  REAL  (G/LB-FORCE) 


Figure  3.  The  Effect  of  Error  on  Three-Dimensional 
Model;  Driving  Point  Real  Acceleration 
Response. 


1* 
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FREE-BODY  ACCELERATION  MOBILITY,  IMAGINARY  (G/LB-FORCE) 


FREQUENCY  -  HERTZ 

Figure  4.  The  Effect  of  Error  on  Three-Dimensional 

Model;  Driving  Point  Imaginary  Acceleration 
Response. 
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ACCELERATION  MOBILITY,  REAL  (G/LB-FORCE) 


Figure  5.  "Beam-Type"  Model;  Free-Body  Real 
Acceleration  Response. 
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ACCELERATION  MOBILITY,  IMAGINARY  (G/LB-FORCE) 


Figure  6.  "Beam-Type"  Model;  Free-Bodv  Imaginary 
Acceleration  Response. 
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are  effectively  invariant  with  error  and  the  method  can  be 
applied,  with  confidence,  to  yield  the  free-body  character¬ 
istics  of  a  structure. 

The  results  of  the  investigation  to  determine  the  influence 
of  measurement  error  on  the  methodology  are  shown  in  Figures 

7,  8,  9,  and  10.  Figure  7  is  a  histogram  obtained  by  con¬ 
sidering  the  real  acceleration  mobility  of  a  specific  element 
on  the  constrained  20-degree-of-freedom  "beam-type"  model 
over  the  frequency  range  of  interest.  The  percentage  error 
was  based  on  a  comparison  of  the  "exact"  or  zero  error  simu¬ 
lation  and  the  conditions  wherein  measurement  error  was 
applied  to  acceleration  mobility  data  and  reaction  force 
data.  The  solid  line  represents  the  theoretical  cumulative 
frequency  polygon  and  yields  the  probability  of  obtaining  an 
error  deviation  below  a  specified  value.  Figure  8  shows 
similar  data  for  the  imaginary  component  of  constrained 
acceleration  mobility.  Corresponding  data  for  the  free-body 
condition  is  given  in  Figures  9  and  10.  Examining  Figures  7, 

8,  9,  and  10  reveals  an  increase  in  error  bandwidth  accom¬ 
panying  removal  of  the  system  constraints.  This  is  a 
consequence  of  polluting  the  reaction  forces  with  measurement 
forces  and  the  subsequent  interaction  with  the  acceleration 
mobility  measurements  also  subjected  to  error.  However,  the 
errors  are  not  prohibitively  large  and  are  within  an  accept¬ 
able  range. 
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PERCENT  ERROR 


Figure  7.  Error  Distribution  lor  Constrained 
Element  Using  Real  Acceleration 
Mobility. 
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PERCENT  ERROR 

Figure  8.  Error  Distribution  for  Constrained 

Element  Using  Imaginary  Acceleration 
Mobility. 
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Figure  9.  Error  Distribution  for  "Beam-n'ype "  Model  Free  Element 
Using  Real  Acceleration  Mobility  Data. 


o 
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PERCENT  ERROR 

10.  Error  Distribution  for  "Beam-Type"  Model  Free  Element 
Using  Imaginary  Acceleration  Mobility  Data. 


CONCLUSIONS 


1.  The  free-body  dynamic  response  of  a  vehicle  can  be  deter¬ 
mined  from  experimental  data  taken  on  a  constrained 
structure. 

2.  The  methodology  is  insensitive  to  measurement  error 
using  simulated  test  data  subjected  to  errors  that  are 
within  the  state  of  the  measurement  art. 

3.  The  method  can  eliminate  the  need  for  soft  suspension 
systems  currently  used  in  helicopter  test  techniques. 

4.  By  eliminating  the  effect  of  the  uncertain  character¬ 
istics  of  the  supporting  structure,  the  accuracy  and 
validity  of  all  dynamic  test  results  will  be  improved. 

5.  The  method  is  inherently  flexible  and  can  be  implemented 
using  displacement,  velocity  or  acceleration  data, 
either  linear  or  angular,  in  two  or  three  directions 

at  one  geometrical  point. 

6.  The  theory  was  found  to  be  experimentally  practical 
and  numerically  sound. 
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APPENDIX 
COMPUTER  PROGRAM 


A  digital  computer  program  was  generated  to  implement  the 
theory  presented  in  this  report.  The  program  was  written 
for  the  IBM  360/40  operating  system  using  FORTRAN  IV 
language.  A  flow  chart  delineating  the  program  logic  is 
shown  in  Figure  11.  A  description  of  the  input  cards  and 
a  program  source  listing  are  included  in  this  appendix. 
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Figure  11  -  Continued. 
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DESCRIPTION  OF  INPUT  CARDS 


Not-*:  All  integer  variables  must  be  right  justified  with 

no  decimal  point. 

Tape,  Card  Reader  and  Printer  Assignments 
1  Card  Reader 
3  Printer 

10  Tape  assignment.  Contains  mobility  data,  with  zero 
error,  for  specified  frequencies. 

All  input  data  must  be  in  the  following  units: 

2 

Mass  -  lb-sec  /in. 

Stiffness  -  lb/in. 

Frequencies  -  Hz 
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PROGRAM  FRESIM 

FREE-BODY  TEST  METHOD  SIMULATION 


Card  1 

Card  2 

Card  3 

Card  4 

Card  5 


Columns  1 

IC 

Program  Control 

2-80 

HEADN 

Heading 

1-5 

NS 

No.  of  stations  with  applied 
force. 

5-10 

INDS 

The  indices  of  the  stations  at 
which  the  forces  are  applied, 
including  all  constraint 
stations  (FORMAT  15,  7110, 
Maximum  of  8) . 

1-80 

SPR 

Spring  rates  at  the  above 
stations,  if  at  constraints, 
otherwise  zero  (FORMAT  8F10.0, 
Maximum  of  8). 

1-80 

FAMP 

Amplitude  of  applied  forces  at 
the  above  stations  (FORMAT 
8F10.0,  Maximum  of  8). 

1-10 

G 

Structural  damping  coefficient 

11-20 

PCT 

Random  error  applied  to  mobility 
amplitude,  uniform  between  - 
and  +  PCT*  element  amplitude. 

21-30 

PCTB 

Bias  error  applied  to  mobility 
amplitude  PCTB*  element 
amplitude. 

31-40 

PHE 

Random  error  in  degrees  applied 
to  mobility  phase  angle.  Uni¬ 
form  between  -  PHE  and  +  PHE. 

41-50 

FPCT 

Random  error  applied  to  force 
amplitude,  uniform  between  - 
and  +  FPCT*  Force  amplitude. 

51-60 

FPCTB 

Bias  error  applied  to  force 
amplitude  FPCTB*  Force 
amplitude. 

61-70 

FPHE 

Random  error  in  degrees  applied 
to  force  phase  angle.  Uniform 
between  -  PHE  and  +  PHE. 

71-80 

IZ 

Random  number  seed. 

if 

w 


/ 
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Card  6  Column  1-5  NFR  Number  of  frequencies  to  be 

used  from  Tape  10  (£100;  if  0 
use  all  frequencies  on  Tape  10). 

6-10  IP1  =  0  Print  constrained  dis¬ 
placement  mobility;  force 
matrix;  free  displacement 
mobility. 

=  1  Print  force  matrix;  free 
displacement  mobility. 

=  2  Print  free  displacement 
mobility. 


Card  7  1-80  INDX  Indices  of  frequencies  to  be 

used  from  Tape  10.  Indices 
must  be  in  ascending  order. 
Five  columns  per  value,  16 
values  per  card  (FORMAT  1615). 
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FREE  BODY  TEST  METHOD  SIMULATION 


USES  MOBILITY  DATA  T-hJM  'XACT' 


INTEGER  HEADNT20I  , HEAO ( 20 } ,HT ( 71  ,1«US(8> ,  !NOX(  100} 

REAL  HZ  (100), YR( 20 .21I.Y! (20t2il  ,iPM HI , F AMP( 8 1 .Y DR ( 20 , 8 » , 

A  YDI(20.8),CYDR(8,9)  ,CYOI  (8,9|  ,MNVK(8.9)  ,F  INVI I  8,  9)  , 

8  FYR(20,8)«FVl(20«8) 

INPUT  HEADING.  XAuf  l APE  HEADING 
100  READ  (1.1101  IC.HEAON 
110  FORMAT  I  ll.A3.19A4) 

IF  (IC.NE.OI  REMIND  10 
IF  (IC.GT.ll  CALL  EXIT 
READ  110)  HT.HEAO.NF.NO. (HZ II  I  , 1"1 , ftF I 
WRITE  (3.120)  H£ADN,HT,HEA0,ND,(HZ(1),I«1,NF) 

120  FORMAT  ( •  1 • /TIO, 17 ( •  FREE  • ) // T2S, Ai, 19A4//T10, *T APE  HEADING'/ 

A  T25.7AA/T25.A3.19AA/T25.I2.'  otoHtfcS  Of  FREEDOM'  /T25, • FREOUENCI  ES 
B(HZ)  ON  TAPE'//(T10.10F10.2I) 

INPUT  SUPPORT  SPRINGS  AND  APPLIED  FORCES 
READ  (1.130)  NS.  1  NOS. S PR .F AMP 
130  FORMAT  1 2 1 5, 71 10/( 8F 10.01 1 

WRITE  (3.1401  1 !NOS(l ),1-1,NS) 

WRITE  13,180)  (  SPRII I  .I-l.NS) 

WRITE  13.1601  ( FAMPI I ) . I • 1 ,NS ) 

140  FORMAT  I /TIO,  'CONSTRAINTS  AND  fUHCtS'//T10,  •  ST  ATI  ON  '.5X.8I10) 
190  FORMAT  ( TIO , • SPR I NG* , 10X ,8F 10.01 
160  FORMAT  ( T 10, • APPL I EO  FORCE  '.•FlO.01 
INPUT  ERRORS,  FREwdcNclES 
REAO  (1,170)  G,  PC T. PCTB , PHE ,F PuT .FPLT8.FPHE, I Z 
170  FORMAT  (  7F10.0.I10I 

WRITE  (3.1801  G 

180  FORMAT  (/TIO, 'STRUCTURAL  OAMP I**' .F6.4I 
IX>IZ*2+1 

WRITE  (3,190)  PCT, PCTB.PHE, IZ ,FPtT .FPCTB, FPHE 
190  FORMAT  (/TIO, 'MAX  RANO  ERROR  BIAS  ERROR  -'.F6.3,'  OF 

ALEMENTS,  MAX  RAND  PHASE  ERROR  «',F5.2,»  DEG.  SEED  -',110/ 

B  TIO, 'FORCE  ERRORS' .4X.F6.3 • 1 9X .F«.3,37X,F9.2/ ) 

REAO  (1,200)  NFR.IPl 
200  FORMAT  ( 161 51 

IF  ( NFR.GT.O)  GO  TO  230 
DO  210  I»1,NF 
210  INDX(II«I 
NFR»NF 

WRITE  (3,220)  NF 

220  FORMAT  (/TlO.'ALL  ',13,'  FREQUtNU  cb  USED* I 
GO  TO  250 

230  REAO  (1  ,2001  ( INOX ( I ) . !■ 1  ,NFR I 

WRITE  (  3,2401  (HZ( I NOXII) 1 , I -  1 ,NFHJ 
240  FORMAT  ( /TIO, ' FREOUENC IES  USE  O' //( f 1U. 10F 12.4) ) 

250  INFR.l 

START  MAIN  LOOP 

DO  410  L'l.NF 


FRE 

1 

FRE 

2 

FRE 

3 

1  FRE 

4 

FRE 

5 

FRE 

6 

FRE 

7 

FRE 

8 

FRE 

9 

FRE 

10 

FRE 

11 

FRF 

12 

FRE 

13 

FRE 

14 

FRE 

15 

FRE 

16 

FRE 

17 

FRE 

18 

FRE 

IS 

;  FRE 

20 

FRE 

21 

FRE 

22 

FRE 

23 

FRE 

24 

FRE 

25 

FRE 

26 

FRE 

27 

FRE 

28 

FRE 

29 

FRE 

30 

FRE 

31 

FRE 

32 

FRE 

33 

fRC 

34 

FRE 

35 

FRE 

36 

FRE 

37 

FRE 

38 

FRE 

39 

FRE 

40 

FRE 

41 

FRE 

42 

FRE 

43 

IFRE 

44 

1FRE 

45 

FRE 

46 

FRE 

47 

FRE 

48 

FRE 

49 

FRE 

50 

FRE 

51 

FRE 

52 

FRE 

53 

FRE 

54 

IFRE 

55 

32 


r>  o 


c 


c 


c 


c 


REAL  1101  FREw, ( ( YRI1 . Jl  .VI  (I  ,J| ,1-l.NU) , J-l.NOI 

1FRE 

56 

IFIL.NE.INOXI  INFRM  GO  TO  410 

1FRE 

57 

ELIMINATE  COLUMNS  Ai«u  CONVERT  TO  ACCL  MOB 

1FRE 

58 

0MR-H2I Ll*6.2832 

IFRE 

59 

OMRS«OMR*OMR 

1FRE 

60 

00  260  l-WNO 

2FRE 

61 

00  260  J-l.NS 

3FRE 

62 

VORII.JI  — VII  1  1 1  NOS  1  J  1  )  *OMR 

3FRE 

63 

260 

VOIII.JI-  VRI 1 « 1  NO  SI J) 1  * OMR 

3FRE 

64 

FORM  CONSTRAINT  DliP  MOB  MATRIX  ANC  ADO  ERROR 

IFRE 

65 

00  270  1-1, NS 

2FRE 

66 

DO  270  J-l.MS 

3FRE 

67 

CVORI  I,  Jl— VORI  INOSII  1  .JI/OMRS 

3FRE 

68 

270 

CVOIII,  Jl— VOIIINOSII  I,  JI/OMRS 

3FRE 

69 

CALL  ERR2  1 VOR , V01 , PC T ,PCTB  ,PHc ,iw,NS ,  IX ,  201 

IFRE 

70 

IF  1  IP1.GT.0I  GO  TO  300 

IFRE 

71 

WRITE  1 3,2801  H2ILI 

IFRE 

72 

280 

FORMAT  l'l',T20, 'SIMULATED  MEA,uittJ  CUNSTRAINEO  ACCELERATION  H0BIL1FRE 

73 

4 1  TV «  REAL  PART  F  m  '.F10.2,*  n2'/l 

IFRE 

74 

CALL  M0UT3  IVOR, NO, NS, 201 

IFRE 

75 

WRITE  (3,290) 

IFRE 

76 

290 

FORMAT  1  /T50, • IMAGINARY  PAKi'/J 

IFRE 

77 

CALL  M0LT3  I VOI ,ND ,NS, 20 1 

IFRE 

78 

ACCELERATIONS  OUE  Tu  FuKLtS  AND  FORCES  AT  CONSTRAINTS 

IFRE 

79 

WITH  ERmJKi  UN  TOTAL  FORCES 

IFRE 

80 

300 

DO  320  J-l.NS 

2FRE 

81 

F-FANPI J| 

2FRE 

82 

00  310  1-1, NO 

3FRE 

83 

YURI  1 , J  I  -YOR 1  I, JI*F 

3FRE 

84 

310 

VOIII.JI- VO III, J)*F 

3FRE 

85 

DO  320  I-l.NS 

3FRE 

86 

Y  YYR  » I—  CYCR I  I  ,  J 1  ♦G*CVOl  1 1  ,  J 1 1  •»•*  jMitl  I  ) 

3FRE 

87 

YYYI-I-CVDI 1 I ,J»-G •CVORI I , J 1  I *r*»P*l  1 1 

3FRE 

88 

CVORI I, JI-YYVR 

3FRE 

89 

320 

CYOHI,  JI-YYVI 

3FRE 

90 

DO  330  I-l.NS 

2FRE 

91 

330 

C YOR 1 1, 1 1 -CVORI 1,1 1 »F  AMP (II 

2FRE 

92 

CALL  ERR2  ICYCR.CVDI.FPCT  ,FPCTo,FPrtt  ,NS,NS,  I X ,  B I 

IFRE 

93 

IF  1  IP1  .GT.ll  GO  TO  360 

IFRE 

94 

WRITE  (3,3401  H2ILI 

IFRE 

95 

340 

FORMAT  1 'l',T20,'SIMULATE0  MEAiJutu  FREE  BODY  FORCE  MATRIX,  REAL 

P  IFRE 

96 

AART  F  *  ' ,F 10. 2  ,  '  H2*/l 

IFRE 

97 

CALL  M0UT3  (C  YOR , NS  , NS, 8  I 

IFRE 

98 

WRITE  1 3,350) 

IFRE 

99 

350 

FORMAT  I  //T50, • IMAGINARY  PrtKl'/l 

IFRE 

100 

CALL  MQUT3  ICVOI , NS, NS, 8 1 

IFRE 

101 

INVERT  rum. t  MATRIX  DIM  8X9 

IFRE 

102 

360 

CALL  CINV2  ICYOR.CVOI  ,NS,FINVR,Fti«Vi  ,8) 

IFRE 

103 

IF  I  IP1  .GT.ll  GO  TO  380 

IFRE 

104 

WRITE  I  3,3701 

IFRE 

105 

370 

FORMAT  I//T40, 'INVERSE  OF  FORCc  rlMlhU,  REAL  PART'/) 

IFRE 

106 

CALI  MQUT 3  ( F I NVR, NS, NS, 8 1 

IFRE 

107 

WRITE  I  3,  3501 

IFRE 

108 

CALL  M0UT3  IF INVI , NS, NS, 81 

IFRE 

109 

FREE  BOOV  MOBILITY  FY  «  YU  •  F  1 NV 

IFRE 

110 
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380  CALL  CMMPY  I  YOR  .  YD  I  ,  F  I  NVR  t  F  I N  V*  ,  itU  ,  NS ,  NS  .  F  YR ,  F  Y  1 , 20 , 8 , 20 »  IFRE  111 

WRITE  (3,3901  H2ILI,(INDS(II,I*i,rol  1FRE  112 

390  FORMAT  ( •  1* , T20 , ‘FREE  600Y  ACC  Moo  MATRIX,  REAL  PART,  F  -SF10.2  1FRE  113 

A  H2'/TS, 101121  1FRE  114 

CALL  M0UT3  IF YR , NO, NS, 201  1FRE  118 

WRITE  13,2901  1FRE  116 

CALL  M0UT3  (F Yl , NO, NS, 201  1FRE  117 

400  INFR» INFR* 1  1FRE  118 

IF  I  INFR.GT .NFR I  GO  TO  100  1FRE  119 

410  CONTINUE  1FRE  120 

REWIND  12  FRE  121 

GO  TO  100  FRE  122 

END  FRE  12? 
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SUBROUTINE  CHUPV  I  A  ,  B  ,  C  ,  0  Nl ,  , ,« j  ,  t ,  F  •  NX  *  ,  MIC  ,  NR  E  I 

C 

COMPLEX  MATRIX  MULT  OBJEcI  1 1  He  DIMENSIONS 
fc  ♦  I«F  » I  A  ♦  I *B) * ( C  ♦  I  *D|  1  «  SQRI-ll 

NR  A  IS  NO  OF  ROMS  IN  DI MENS* Ju  JF  A.B 

NRC  IS  NU  CF  ROMS  IN  OIMENS4U1.  JF  CiD 

NR E  IS  NO  OF  ROMS  IN  DIMENStu*  jF  t,F 

A,B  ARE  Nl  X  N2  C.O  ARE  >«e  A  NJ  E,F  /.RE  N1  X  N3 

REAL  A)  NRA, n  tBINRA.ll ,C I NRC  .  1 1  ,  ji NKu , 1)  ,  E I  NRE  .  11  tFINRE.l > , 

A  G (20, 20) 

CALL  MMPY2  I A,C,N1 ,  N2  ,  N3  ,  E  ,  NR  A  ,  ,kl  ,  NKc.  > 

CALL  MMPY2  IB.D.Nl , N2 , N3 . G > NR  A , unc f eu) 

DU  ICO  I  -  1 . Nl 
DO  ICO  J«1.N3 
100  Ell . J)*E( I , JI-GI I , Jl 

CALL  MMPY2  IA,0,Nl,N2,N3,F,NRA,,<(u.,KRtl 
CALL  MMPY2  IB»C»N1«N2,N3»G«NRA,  t4.0  I 
00  110  I ■ l *  N 1 
00  110  J* 1 »N3 
110  FII.J)*F< I . J I «G I  1 . J ) 

RETURN 
END 


CMY 

l 

CM  ’ 

2 

CM, 

3 

CMY 

A 

CMY 

X 

CMY 

6 

CMY 

7 

CMY 

8 

CMY 

9 

CMY 

10 

CMY 

11 

CMY 

12 

CMY 

13 

1C  MY 

1  A 

2CMY 

15 

2CMY 

16 

CMY 

17 

CMY 

lfl 

1CMY 

19 

2CMY 

20 

2CMY 

21 

CMY 

22 

CMY 

23 
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SUBROUTINE  CINV2  ( * , B,N,C ,D,NR| 

CIN 

1 

c 

CIN 

2 

c 

COMPLEX  INVERSE  C  ♦  1*0  -  1.4V  Jl  1  »  l*B 

I  -  SQR  1-1) 

CIN 

3 

c 

OBJECT  TIME  DIMENSIONS  NR  la  .41*  uF  ROWS 

IN  OIM.  OF  A , B ,C , D 

CIN 

A 

c 

THEY  MUST  BE  DIMENSIONED  Mlfrl  caIHA  COLUMN 

IN  CALLING  PROG. 

CIN 

5 

c 

CIN 

6 

c 

A  IS  ASSUMED  NON-SINGULAR  N 

MUST  BE  .LE.  20 

CIN 

7 

c 

USES  INVRS2.  MMPY2 

CIN 

B 

c 

CIN 

9 

REAL  AINRtD.BINRtll.CINRtll.OlNMf  U.fc(20«2i: 

CIN 

to 

CALL  INVRS2  1C, N, A, NR) 

CIN 

11 

CALL  MMPY2  IC,B,N,N,N,E,NR»NR»2Ul 

CIN 

12 

CALL  MMPY2  (B,E,N,N,N,D,NR,20,i4KI 

CIN 

13 

DO  100  I«1,N 

1C  IN 

1 A 

DO  100  J-1,N 

2C  IN 

15 

El  1 ,J)«-EII,J) 

2C  IN 

16 

100  Dl I , JI-OI 1 , J)«AII  ,  J) 

2C  IN 

17 

CALL  INVRS2  IC,N,D,NRI 

CIN 

IB 

CALL  MMPY2  IE,C,N,N,N,D,20,NR,.4K) 

CIN 

19 

RETURN 

CIN 

20 

END 

CIN 

21 
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100 


SUBROUTINE  MMPY2  I  A ,B , N1 ,N2 .N3 ,t ,NHA» NRBt NRC) 

MV2 

1 

OBJECT  TIME  DIMENSIONS  NRm  ii 

NU  OF  ROMS  IN  A  IN  DIMENSION 

MY2 

2 

STATEMENT  I K  MAIN  PROG.  ETC  *-UK 

NKb»  NRC. 

MY2 

3 

C  ■  A  •  b 

MY2 

A 

A  INI  X  ,<Ul 

0  IN2  X  N3)  C  INI  X  N3I 

MY2 

! 

MY2 

6 

REAL  AINRA,  ll.BINRB, 11. CINRC.ll 

MY2 

7 

00  ICO  I •  1  «N1 

1MY2 

8 

00  ICO  J«i,N3 

2MY2 

9 

C ( » # JI*U. 

2MY2 

10 

00  ICO  K* 1  * N2 

3MY2 

11 

C(  I  ,  J)*C(  I.JMAU  J) 

3MV2 

12 

RETURN 

MY2 

13 

END 

MY2 

l  A 
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SURR00T  1NE  RANDU  (  1 X ,  I  V •  YFL  1 

RAN 

1 

THIS  SUBROUTINE  la  l-nUH  SSP  VERS.  II 

RAN 

2 

IV»I XP65539 

RAN 

3 

iHiviiooaio.no 

RAN 

A 

100 

I  V» I V*21A7483647*l 

RAN 

*i 

110 

VFL-IV 

RAN 

6 

VFt »VFL *• A656613E-9 

RAN 

7 

RETURN 

RAN 

B 

END 

RAN 

R 
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SUBROUTINE  ERR2  ( A  ,B .  PCT  ,PCTB  Ml,N2,  IX .  NR  I 

ERR 

1 

OBJECT  TIME  DIMENSIONS 

ERR 

2 

EACH  ELEMENT  OF  A  COMPLEX  MATRIX.  A  ♦  l*B.IS  MODIFIED  TO 

ERR 

3 

INCLUDE  A  SMALL  PHASE  cKKUM •  PNE  IOEG).  A  BIAS  ERROR. 

ERR 

4 

PCTB  (RATIO!  ON  AMPLITUDE.  AND  A  UNIFORM  RANDOM  ERROR 

ERR 

5 

HAVING  A  ♦/-  MAXIMUM  OF  FcT  (RATIO!  ON  AMPLITUDE. 

ERR 

6 

THE  PHASE  ERROR  IS  ALSJ  HAi.uUMLY  DISTRIBUTED 

ERR 

7 

TE  NO  SYMMETR IZAT1QN  IS  PERFjhMcJ 

ERR 

B 

ERR 

9 

USES  RANOU 

ERR 

10 

DIMENSION  A(NR«1) .BINR.ll 

ERR 

11 

IFIPCT  1  120.100*120 

ERR 

12 

100 

IF  1 PCTB )  120,110.120 

ERR 

13 

no 

IFIPHEI  120.140.120 

ERR 

14 

120 

P-P ME/57.296 

ERR 

13 

DO  130  I-1.N1 

1ERR 

16 

DO  130  J*l,N2 

2ERR 

17 

CALL  RANOU  (IX.IY.YFL! 

2ERR 

IB 

IX*  1 Y 

2ERR 

19 

E*2. 0*P*( YFL-O. 51 

2ERR 

20 

Al-AI  1. J!-E*B(l,Ji 

2ERR 

21 

B  (  I  ,  J  )■  B(  I  ,  J  )  »E*A(  I  ,  J  ) 

2ERR 

22 

A( I, Jl-Al 

2ERR 

23 

CALL  RANOU  (IX.IY.YFL! 

2ERR 

24 

IX*!  Y 

2ERR 

23 

E*l.0*2.0*PCT*IVFL-0.5|4PCTB 

2ERR 

26 

A ( 1 . JI«AI I, JI*E 

2  ERR 

27 

130 

B( I . J!«B(  I,  J»*E 

2ERR 

26 

140 

RETURN 

ERR 

29 

END 

ERR 

30 

I 

i 
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SUBROUTINE  MOUT3  (A,M,N,NR) 

C  OBJECT  TIME  DIMENSIONS  Nh  Ij  NO  OF  OIMENSIONEO  ROMS 

C  N  SHOULD  NOT  EXCEEO  20 

REAL  A(  NR , l ) 

ID*MINO(N» 101 
WRITE  I  3.1001  I  1.1*1.101 
100  FORMAT  I/T5, 101121 
WRITE  (3.100) 

DO  110  I  * l.M 

110  WRITE  (3,1201  I.Ud.JI.J-l.lD! 

120  FORMAT  (I5,5X,IP10E12.*) 

IF  (ID-N)  130, ISO, 150 
130  WRITE  (3.100)  (1,1*11. N) 

WRITE  (3,100) 

DO  1  AO  1*1, M 

1*0  WRITE  (3,120)  I . ( A ( I • J ) . J*1 1,N) 

ISO  RETURN 
END 


M03 

1 

M03 

2 

M03 

3 

M03 

* 

M03 

5 

M03 

6 

M03 

7 

M03 

a 

1M03 

9 

1M03 

10 

M03 

11 

M03 

12 

M03 

13 

M03 

1* 

1M03 

15 

1M03 

16 

M03 

17 

M03 

18 
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SUBROUTINE  INVRS2  tA.N.B 

,NR) 

INV 

OBJECT  TIME  DIMENSIONS. 

mH KM/4  SHOULD  BE  DIMENSIONED 

INV 

W  1  TH  ONE  MORE  COLUMN 

THAN 

Run. 

NR  SHOULO  BE  .LE.  20 

INV 

A  -  INVERSE  OF  B 

B  UNDI jfUKjtL 

INV 

INV 

01 MENSI ON  A(NR,l),D(20«21),IR0al2i), ICOL  1211, BI NR, l) 

INV 

00  100  I » l »N 

1  INV 

00  ICO  J»1  »N 

2  INV 

100 

At  I , J)*BI l,J) 

2  INV 

M»NM 

INV 

00  110  l  »l,N 

1  INV 

IRON  1 I ) « I 

l  INV 

110 

ICOL  <11*1 

1  INV 

DO  2 60  K«1,N 

1  INV 

AMAXa  A  IK,K) 

1  INV 

00  130  I-K.N 

2  INV 

00  130  J*K,N 

3INV 

IFIABSI  A1I»JI)-ABSIAMAX)I130 

.ICO,  120 

3  INV 

120 

AMAX>  At  I.JI 

3IN7 

1C  *  I 

3  INV 

JC-J 

3INV 

130 

CONTINUE 

3  INV 

KI«IC0L(KI 

1  INV 

ICOL  t  K 1  a ICOLI IC1 

1  INV 

ICOL  t  IC  )  —  K I 

1  INV 

K I  - ) ROM  IK) 

1  INV 

IROUIK) a I  ROW 1 JC I 

l  INV 

IKOWI JC  >»KI 

1  INV 

l  F  (  A  MAX  1  160.140,160 

l  INV 

140 

WRITE  (3,150) 

1  INV 

1  30 

FORMAT!  •  SOLUTION  OF  EXISTING 

NAIttiX  NOT 

POSSIBLE') 

1  INV 

GO  TO  330 

1  INV 

160 

CO  170  J>1,N 

2  INV 

E*AIK,J) 

2  INV 

AIK, JlaAt IC , J 1 

2  INV 

170 

At IC.JI -E 

2  INV 

DO  1E0  I ■ 1 , N 

2  INV 

EaAt  l,K> 

2  INV 

A(I,K)aA(I,JCI 

2INV 

180 

All, JC) «E 

2  INV 

DO  210  I  a  1 , N 

2  INV 

IF(I-K)  200,190,200 

2  INV 

190 

At  I ,M|«1. 

2  INV 

GO  TO  210 

2  INV 

200 

At  I ,N)»0. 

2  INV 

210 

CONTINUE 

2  INV 

PVT=A(K ,KI 

1  INV 

DO  220  J ■ l , M 

2  INV 

220 

A(K, J)«A(K,J)/PVT 

2  INV 

00  250  1*1, N 

2  INV 

IFt I-K) 23C.25C.230 

2  INV 

230 

AMULT’A  < 1 , K 1 

2  INV 

00  240  J* 1 , M 

3  INV 

240 

A(l,J)aA(I,J)-AMULT*A(K, 

J) 

3  INV 

250 

CONTINUE 

2  INV 

V 
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00  260  I ■  1, N 
260  Af I,K»*A(l,M) 

00  290  I«1,N 
00  270  L«1,N 

IF  ( 1ROWII1-U  270,280, 270 
270  CONTINUE 
280  00  290  J*  l ,  N 
290  D(L,J(»A( I , J) 

00  320  J»1,N 
00  300  L-I,N 

IF!  ICOL (JI-LI  300,310,300 
300  CONTINUE 
310  00  320  1*1,  N 
320  AI!,l)-OII,JI 
330  RETURN 
END 


2INV 

56 

2  INV 

57 

1INV 

58 

2  INV 

59 

2  INV 

60 

2  INV 

61 

2  INV 

62 

2  INV 

63 

1  INV 

64 

2  INV 

65 

2  INV 

66 

2  INV 

67 

2  INV 

68 

2  INV 

69 

INV 

70 

INV 

71 
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